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ABSTRACT 

Using observations from the Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE) 
survey of the Magellanic Glouds, we have found thirty five evolved stars and stellar end products 
that are bright in the far-infrared. These twenty eight (LMG) and seven (SMC) sources were se¬ 
lected from the 529 evolved star candidates in the HERITAGE far-infrared point source catalogs. 
Our source identification method is based on spectral confirmation, spectral energy distribution 
characteristics, careful examination of the multiwavelength images and includes constraints on 
the luminosity, resulting in a thoroughly vetted list of evolved stars. These sources span a wide 
range in luminosity and hence initial mass. We found thirteen low- to intermediate mass evolved 
stars, including asymptotic giant branch (AGB) stars, post-AGB stars, planetary nebulae and 
a symbiotic star. We also identify ten high mass stars, including four of the fifteen known B[e] 
stars in the Magellanic Clouds, three extreme red supergiants which are highly enshrouded by 
dust, a Luminous Blue Variable, a Wolf-Rayet star and two supernova remnants. Further, we 
report the detection of nine probable evolved objects which were previously undescribed in the 
literature. These sources are likely to be among the dustiest evolved objects in the Magellanic 
Clouds. The Herschel emission may either be due to dust produced by the evolved star or it may 
arise from swept-up ISM material. 

Subject headings: stars: late-type - infrared: stars - Magellanic Clouds - circumstellar matter - stars: 
mass-loss - - submillimeter: stars 
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1. Introduction 

Evolved stars are important drivers of the 
chemical enrichment of the interstellar medium 
(ISM) of galaxies. In our own Galaxy dust enrich¬ 
ment of the ISM is dominated by the dusty winds 
from low- to intermediate mass stars (1-8 Mq) 
during the as ymptot i c gian t branch (AGB) phase 
of evolution ( Gehrz 1989ll . At lower metallici- 


ties the relative contribution from high mass stars 
(M > 8Mq) becomes important, and a significant 
period of dust production from these high-mass 
sources, via their stellar winds or directly in the 
supernova (SN) ejecta, is often invoked to ex- 


in high-redshift galaxies and 

auasars (IDwek et al. 

2007; Michalowski et al. 20It 

t Beelen et al. 

20061 

Bertoldi & Gox 

2002|1. 


Dust production at low-metallicities is only 
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just beginning to be understood, due primar¬ 
ily to recent galaxy-wide surveys of nearby low 
metallicity systems such as the Magellanic Clouds 
(MC). The Large and Small Magellanic Clouds 
(LMC and SMC, respe ctively), with metallici- 
ties of 0.5 and 0.2 Zp, (Russell fc Dooital Il992 : 
iTchernvshvov et all 2015 1. are ideally suited to 
such inv estigations, due to their p r oximity (50 and 
60 k pc: Ngeow fc Kanbur 20081 : ISzewczvk et al 


200911 ■ favorable viewing ang le (35° for the LMC; 
van der Marel fc CionillioOlIl and low average red- 
deni ng along the line-of- sight {Eb-v ^ 0.08 and 


0.04; Schlegel et al. 1998ll . 

The resolved stellar populations of the Magel¬ 
lanic Clouds have been well studied over a wide 
wavelength range from the optical to the far- 
infrared (e.g. Zaritskv et al. llOOTtl iHarris fc Zaritsk"^ 
|2004 I Blum et ~ 2006 : Bover et al. 201lh , as 
these surveys cover the full spatial extent of each 
galaxy, a complete picture of the mass-losing stel¬ 
lar populations in the Magellanic Clouds can been 
obtained. The dust-injection rate from low-mass 
stars in the Magellanic Cloud s has been esti¬ 
mated bv Srinivasan et al. ( 200^ Matsuuraet_^ 


( 20091) : Riebel et al. (2012); Bover et al.l ( 2012 ). 


using mid-lR fluxes which are sensitive to warm 
100 K) dust. To compile a complete inventory 
of dust producers and fully characterise the stel¬ 
lar dust return to the ISM, the far-infrared (IR; 
50-250 fim) and sub-millimeter (submm; (A ^250 
^m) emission from cold dust (T < 50 K) envelopes 
associated with stars in an evolved stage of stellar 
evolution must also be constrained. 

Using data from the Herschel Inventory of the 
Agents of Galaxy Evolution in the Magellanic 
Clouds Open Time Key Project (HER ITAGE; 
iMeixner et al.l l2013l) . Seale et Si ( 20141) created 


a Herschel band-matched source catalog (BMC) 
for both the LMC and SMC by combining the 
individual HERITAGE catalogs of point sources 
at 100, 160, 250, 350, and 500 /rm. They clas¬ 
sify the sources in the BMG and find that the 
brightest HERITAGE sources are predominately 
young stellar objects (YSOs) while the faintest are 
background galaxies. Intriguingly, they provision¬ 
ally identify a smaller population (^500) of very 
dusty stars, in the later stages of stellar evolu¬ 
tion (e.g. extreme AGB stars, post-AGB (PAGB) 
objects, planetary nebulae (PNe) and supernova 
remnants (SNR)), although the majority of these 


sources have also been identified as YSOs in other 
studies. 

In the Milky Way, the Hersch el Key Program 
MES S (Mass loss of Evolved StarS: lGroenewegen et al, 
201l[) . has observed a representative sample of 
low- and intermediate-mass and high-mass post- 
main sequence objects, across a range of mass- 
loss rates. They find that extended structure 
in the far-IR is comm only seen around AGB 
stars and their remnants (Kerschbaum et al 1I2H^ 

Cox et al. 2012 : iMaercker et al. 2012 ). However, 
the presence of a significant number of evolved 
stars detected in the HERITAGE survey is un¬ 
expected, as it was predicted that the major¬ 
ity of these sources in the Magellanic Clouds 
would fall below th e Herschel sensitivity limit 


( Meixner et al. 20101) . Only the most luminous. 


heavily dust enshrouded evolved stars were ex¬ 
pected to have a detectable submm emission. 
The cold circumstellar envelopes of two mas¬ 
sive evolved stars (R71 and IRAS 05280-6910) 
i n the LMC have been detected with Herschel 
( Bover et al. 20101 ). However, I van Loon et al.l 


(|2010bl ) attribute this emission to swept-up ISM 
dust rather than from dust produced by the stars 
themselves. 

The goal of this paper is to investigate the 
^500 evolved star candidates with a submm ex¬ 
cess in the Magellanic Clouds identifie d using data 


from the Herschel HERITAGE-BMC (Seale et al 


20141) . The paper is organised as follows: the 
observational data is summarized in Section 
In Section |3] we describe the selection of post- 
main-sequence stars from the HERITAGE-BMC; 
in Section 0] we present the results of this search 
and discuss the different post-main-sequence ob¬ 
ject classes contained in the catalog. Einally, Sec¬ 
tion [5] concludes this paper. 

2. Photometric Data 

Ear-infrared (far-lR) and submm observations 
of the Magellanic Clouds were taken as part 


of th e HERITAGE Key Project (jMeixner et al 


26131) using the Photode tector Array Camera 
and Spectrometer (PACS; IPoglitsch et al.l 20161) 
and the Spectr al and Photometric Imaging Re¬ 
ceiver (SPIRE; iGriffin et al.l 2010) on the Her¬ 
schel Space Observatory ( PilbrattF 2010l) . Photo¬ 
metric catalogs of point sources were created us- 
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ing STARFINDER ( Diolaiti et al. 20001) which effec¬ 
tively removes diffuse emission to extract a point 
source flux. A full description of the HERITAGE 
observat ions and data reduc tion details can be 
found in Meixner et al. ( 20131) . 

The HERITAGE survey detected over 35000 
unique sources in the EMC and 7500 in the SMG 
at 100, 160, 250, 350, and 500 /rm. T hese sources 
were positionally cross-matched by ISeale et ^ 


( 20141) to produce a HERITAGE Band-Matched 
catalog (BMG) for the LMG and SMC. This cata¬ 
log was further matched to the Spitzer IRAC and 
MIPS catalogs from the Surveying the Agents 
of Galaxy Evolution (SAGE) Spitzer legacy sur- 
yeys of the Large and Small Magellan ic Clouds 
( Meixner et al.l[2C)06HGordon et al.feoilh . Eurther 
details rega r ding the band-matching are given in 
Seale et al.1 (|2014l) . The final HERITAGE-BMC 


covers the 3.6 to 500 fim range and contains over 
42000 unique sources; these represent the dustiest 
populations of sources in the Magellanic Clouds. 

In order to comprehensively classify the 500 
sources in the HERITAGE photometric cat¬ 
alogs, and fully explore the dust production 
from point sources in the Magellanic Clouds, 
we have cross-matched the HERITAGE-BMCs 
with optical, near-IR and mid-IR point source 
catalogs, to provide photometry over a broader 
wavelength range (0.36 - 500 /rm). The opti¬ 
cal UBVI and near-IR JHKg photometry comes 
from the M agellanic Clouds Ph otometric Sur¬ 
vey (MCPS; IZaritskv et al.l [l9^ : the Infrared 
Survey Eacility (IRSE) Magella nic Clouds Point 


Source Survey (jKato et al. 


2007 ), the Two Micron 


All-Sky Survey (2MASS; ISkrutskie et al.l I 2 OO 6 I) . 
and the 2M ASS 6X Deep Point So urce Catalog 
(6X2MASS; ICutri fc 2MASS TeamI l2004 . We 
also obtained supplementary mid-IR d ata from 


the AKARI LMC poin t source catalog (jlta et al 


l2008l : lKato et al.l2007 ) and the WISE all-sky cat 
alog (|Wright et al. 20I0l) . 

We considered sources from the various catalogs 
to be a match if the source is within 5" of the IRAC 
3.6 fim position in the BMC. In cases where there 
were multiple matches between catalogs we only 
considered the closest match. 

To complement the photometric data we searched 
the spectroscopic data from the Sp itzer SAGE- 


Spec survey of the Magellani c Clouds (|KemDer et al. 


2 OIOI : yan Loon et al. 2010b ) for all spectra taken 


with the InfraRed Spectrograph (IRS; 5.2-38 /rm; 
Houck et al. 2004) and MIPS spectra (MIPS-SED; 
52-97 fim) for matches to the HERITAGE-BMCs. 
We found 185 matches. These sources have 
been s pectroscopica l ly cla ssified by IWoods et al.l 
(I 2011 I) : I Ruffle et ah ( 2015 ) and Woods et al. (in 
prep). 


3. Refining the identification of post main 
sequence stars 


To identify the evolyed stellar population in the 
Magellanic Clouds with a Herschel detection, we 
compiled a preliminary candidate list from sources 
which were identified as dusty objects in the late 
stages of stellar eyolution including extreme and 
post-AGB stars, planetary nebulae, and super- 
nova remnants. T hese sources were categorized by 
ISeale etliL ( 20141) by matching the HERITAGE 
catalogs to other classification schemes and lists 
of preyiously identified objects in the literature. 

The HERITAGE BMCs contain 443 LMC and 
51 SMC evolyed star ca ndidates and 45 L MC plus 
9 SMC PN candidates ( Seale et al. 2014 ). In to¬ 
tal 455 of these Magellanic Cloud sources have also 
been assigned an alternate classification, making it 
difficult to judge the reliability of the source clas¬ 
sifications in the initial list of evolved stars. For 
this reason we have to employ additional selection 
criteria to remove young stellar objects and back¬ 
ground galaxies from the evolved star candidate 
list. 


The following subsections describe our selection 
process, the methodology is generally listed in or¬ 
der of procedure however at some stages the se¬ 
lection criteria were assessed in parallel. For in¬ 
stance, an unambiguous identification can not be 
obtained solely upon the inspection of the source’s 
spectral energy distribution (SED), thus the ob¬ 
ject’s position and morphology was also considered 
simultaneously when assigning a likely classifica¬ 
tion. 


3.1. Spectroscopic confirmation 

The most reliable point-source classifications 
are obtained from spectroscopic observations. 
This alleviates most of the ambiguity that arises 
from photometric classification schemes due to 
overlaps in colour-magnitude space and thus pro¬ 
vides robust source classifications for IR objects. 
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We compare our initial list of evolved star candi¬ 
dates to the SAGE-Spe c spectroscopic surv ey of 


the Magellanic Clouds ( Kemoer et al. 201011 . Of 


the 529 unique sources in our list, 172 were found 
to have a spectroscopic counterpart. 

From IRS spectroscopy we find a total of 
seven evolved stars and four planetary nebulae. 
Other less-common objects were also identified; 
these include one blue supergiant (BSG), one 
B[e] star, one Luminous Blue Variable (LBV), 
one Wolf-Rayet (WR) star, one symbiotic star 
and two sources of an unknown phys ical nature 
( Woods et al. 201l[ iRuffle et al.ll2015i Woods et 
al. in prep.). Eight of these sources also have 
far-IR MIPS s pectra covering the 52- 97/rm wave- 


spe __ 

length range (Ivan Loon et al.ll2010bl la[l. The re¬ 
maining 154 objects with an IRS counterpart 
were co nfirmed spectroscop i cally as YSOs or Hii 


regions (IWoods et al.l 12011 : Oliveira et al. 2013 : 


Ruffle et al. 20151) . We remove these sources from 


the list of the evolved star candidates. The re¬ 
sults of the spectral comparisons are summarised 
in Tables HI and m 


From the IRS spectra we have identified 18 
post main sequence stars. Spectroscopic confir¬ 
mation is still needed for the majority 60%) 
of the BMC evolved star candidates. For these 
357 sources, we use several other complementary 
metrics in parallel to find evolved objects, improve 
upon the source identification and to exclude non- 
evolved contaminants from our sample. 


3.2. SED classification 


Photometrically it can be difficult to separate 
AGB stars from YSOs and young PNe from com¬ 
pact Hii regions. The colour select i on cri teria for 
the far-IR sources from iBover et a . J 20n ) and the 


candid ate post-AGB stars from Ivan Aarle et al 
( 201ll) heavily overlap with the same regions as 
the YSO candidates in the colour magnitude dia¬ 
grams (CMD). As such many of these sources have 
multiple classifications in the literature and we ex¬ 
pect a high degree of contamination from YSOs in 
our initial evolved star sample selection. 


To refine our selection further, the nature of 
each object with a far-IR detection was assessed 
by examining its SED from near-IR to far-IR wave¬ 
lengths. FigurelHshows some illustrative examples 
of typical SEDs for the different classes of objects. 


Based on the shape of the SED with the flux plot¬ 
ted in log(AFA) against log(A) and more specif¬ 
ically on the peak of the energy distribution we 
can (tentatively) divide the sources in our list into 
four separate classes: (1) evolved star candidate, 
(2) probable YSO, (3) possible background galaxy, 
(4) a mismatch between the near- and far-IR data 
used in this study. 


3.2.1. Evolved star candidates 


The most extreme AGB stars can have red 
colours similar to YSOs, background galaxies, and 
PNe, and will also display a rising SED between 8 
to 24 /rm. We identify dust enshrouded AGB stars 
based on the peak of the SED; dusty AGB stars 
typically peak in the IRAC bands (<8 /rm), while 
stars with little to no circumstellar dust peak in 
the near-IR. As the AGB star becomes more ob¬ 
scured the peak of the SED shifts to longer wave¬ 
lengths, and the SED may peak between 8 and 
24 /im for the most extreme AGB stars. AGB 
stars detected with Herschel will fall in the lat¬ 
ter category due to the cool dust envelopes which 
are likely providing additional extinction at the 
shorter wavelengths. In all cases the shape of the 
SED for AGB stars resembles that of a blackbody 
with temperature between 100 and 3000 K. 


Red supergiants are more massive and are in 
general more luminous than AGB stars. As they 
exhibit a similar dust chemistry to oxygen-rich 
AGB stars, their SEDs are virtually indistinguish¬ 
able. Ancillary information such as luminosity, 
pulsation period or membership of a young cluster 
is generally required to separate these classes (see 
Section 13.41) . 

Post-AGB stars tend to have a double-peaked 
SED, with one peak in the optical/near-IR (due 
to a hot central stellar component) and one in 
the mid-IR (due to a detached shell of circum¬ 
stellar emission). This double-peaked morphology 
enables one to readily distinguish post-AGB stars 
from SEDs with a stellar component which has a 
large infrared excess. 


The SEDs of other more massive evolved stars 
such as LBVs, WR stars and BSGs tend to be 
double-peaked, showing a stellar photosphere and 
a strong far-infrared excess. The classifications 
for these sources are taken f rom the literature 
(e.g. Bonanos et al.llioO^ 201(ll) . and are not solely 
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based on the SED. 


3.2.2. YSO candidates 



Wavelength |/<m] 

Fig. 1.— Example SEDs for the different groups 
of objects, including a YSO, a background galaxy, 
an AGB star with a large dust excess, a double- 
peaked SED indicative of a PAGE star, and a PN. 
The shape of the SED for YSOs differs markedly 
from the evolved stars. 


The SEDs of Stage O/I YSOs (early to mid¬ 
stage YSOs) show a monotonic increase from the 
mid-IR to the far-IR, with the peak at ^100 /rm. 
Objects which peak in the far-IR are heavily ob¬ 
scured by cold dust, which absorbs the photo- 
spheric emission from the forming star, and re¬ 
radiates this energy at far-IR wavelengths. As 
the YSO evolves it dissipates its surroundings and 
emission from warm dust begins to dominate the 
SED. For stars in the latest stages of formation, 
the peak emission is blue-ward of the far-IR, and 
these objects can be difficult to disentangle from 
other so urces with a warm c ircumstellar dust con¬ 
tinuum ( Oliveira et al Hiooi . Given the high like¬ 
lihood of contamination from YSOs in our sample, 
additional care needs to be taken when inspect¬ 
ing the shape and peak of the energy distribution 
for these objects. However, objects which have an 
SED peak between 70 a nd 150 /xm are likely YSOs 
(e.g. Sewilo et al 


To break the degeneracy between YSOs in the 
more advanced stage of formation and evolved 
stars which exhibit a significant IR excess, we cre¬ 
ated a mean YSO SED for each evolutionary stage 
based on th e 154 YSOs with a spectroscopic coun¬ 
terpart Isee lWoods et al.lfeoill for each YSO spec¬ 
tral class). These YSOs range from Stage 0/1- 
4 and cover a wide range in luminosity. Using 
the YSO SED templates as a basis, we only re¬ 
tained those objects which peaked at wavelengths 
< 70^m or had a significantly different profile to 
that of a YSO. We find 197 evolved star candi¬ 
dates closely mimic the SEDs of the previously 
identified YSOs. After their removal, our sample 
contains 178 dusty evolved star candidates. 


3.2.3. PNe & background galaxies 

The SEDs of PNe have a very similar appear¬ 
ance to those of YSOs, both in terms of shape and 
in the position of the peak; rising from the mid- 
IR and peaking in the far-IR. Background galaxies 
also tend to have a rising SED in the mid-IR, and 
the peak of their SEDs tends to lie in the Her- 
schel bands. We can not resort to SEDs to distin¬ 
guish between these classes. We further consider 
potential confusion with background galaxies in 
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Section 13.3.21 

3 . 2 . 4 . Mis classification and mismatches 

At low flux levels where a source is only de¬ 
tected in a few Herschel bands, the shape of the 
SED can be difficult to characterize. This can also 
be problematic in crowded fields, such as star¬ 
forming regions where there is a high probabil¬ 
ity of misidentification between photometric cata¬ 
logs, due to the wide range in resolution between 
the surveys. We examine the SEDs for any large 
discontinuities between the photometry from dif¬ 
ferent instruments and point source catalogs, in¬ 
dicative of a likely mismatch. An example of 
a mismatch between the photometric catalogs is 
shown in the bottom panel of Figure [5J The SEDs 
for these sources are most likely a combination 
of a point source with warm circumstellar dust 
and spatially coincident starless dust clump which 
radiated in the far-IR. We eliminate 109 BMC 
evolved star candidates from the sample due to 
a discontinuity in the SED. 

Low flux levels may also be problematic when 
comparing sources in the HERITAGE-BMC to ob¬ 
jects previously identified in the literature, as this 
leads to higher positional uncertainties. This po¬ 
sitional uncertainty in turn results in the potential 
misclassification of the BMC source as an evolved 
star. It is therefore vital that any source identi¬ 
fied in the HERITAGE-BMC as an evolved star 
has a smooth SED with data points that appear 
consistent with each other and are associated with 
a point source. 

3.2.5. SEDs of the final candidate list 

Figures [3] and [7] show the SEDs (as well as the 
available spectral data) for the final candidate list 
of post main-sequence objects in the MCs. In all 
instances a smooth transition is seen between the 
2MASS, Spitzer and Herschel photometry, thus 
the shapes of the SEDs in the IR are well defined. 
At shorter wavelengths it is sometimes unclear if 
the MCPS data is associated with the source due 
to the difference in resolution, the high circumstel¬ 
lar extinction at optical wavelengths and intrinsic 
variability of the source. 

The SEDs show a range of distinguishing fea¬ 
tures due to the presence of warm and cold dust, 
which is dependent on the mass-loss history of the 
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Fig. 2.— An example SED for a source in the 
BMC which is similar in appearance to that of a 
YSO (top) and for a source where there is a mis¬ 
match between the Herschel flux and the photo¬ 
metric data at shorter wavelengths (bottom). 
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source. In all instances the SEDs are consistent 
with post main-sequence stars; the final classifica¬ 
tion is noted in each panel. 


3.3. Image inspection 


Using multiwavelength observations we assess 
the nature of each source in our initial sample. In 
several instances a source may only be identified 
in one or two Herschel bands and may be missing 
counterparts at other wavelengths. Simutanious 
inspection of the images and the SEDs ensures 
that the BMC sources are real and not a com¬ 
posite of several nearby sources which contribute 
to the emission at different wavelengths, resulting 
from a mismatch between the Herschel and Spitzer 
catalogs. 


The PACS images are known to suffer from 
low level artifac ts (in the form of striping and 
cross-hatching) (jMeixner et al.l 120131) which may 
affect the quality and reliability of point sources 
extracted from the original images at 100 and 
160/rm. Thus care needs to be taken when study¬ 
ing compact objects detected in the PACS bands. 

In an effort to obtain high-quality PACS im¬ 
ages at 100 and 160^m, the HERITAGE PACS 
data have been reprocessed with the Tamasi^ 
software (Chanal and Panuzzo, private communi¬ 
cation) which mitigates thermal drifts in the PACS 
imaging using the redundancy and high pass filter¬ 
ing of the observations, effectively removing the 
artificial structures that compromised the origi¬ 
nal images. However, the Tamasis images cannot 
be used to extract photometry because the map¬ 
ping technique and filtering makes the photometry 
potentially systematically inaccurate. Thus, we 
must inspect the refined PACS Tamasis images to 
distinguish between striping and point-sources in 
the HERITAGE BMG. We remove three spurious 
sources that are only detected with PACS and are 
coincident with the striped/hashed regions in the 
original images. 


3.3.1. Evolved stars vs. YSOs 

To determine the likelihood that a given source 
is an evolved star we conduct a careful vi¬ 
sual inspection of the IRAC, MIPS, PACS and 
SPIRE images. Magellanic Cloud Emission Line 


^http://pchanial.github.io/tamasls-pacs/ 
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Fig. 3.— SEDs and spectra for the post-main se¬ 
quence sources in the LMC, ordered in terms of 
RA. Data from the HERITAGE-BMC and SAGE 
catalogs are shown as the dashed-diamond line, 
MGPS data as crosses, WISE data as triangles 
and AKARI data as squares. 
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Fig. 4.— Fig. |3l Continued. 


Fig. 5.— Fig. |3l Continued. 
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Fig. 6.— Fig. |3l Continued. Fig. 7.— As for Fig.[3l but now for the SMC. 
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Survey (MCELS; ISmith fc MCELS TeamI 1998 ) 
Hg images wer e also examined. Using DS9 
( Jove fc Mandeil 1200,81) we simultaneously view 
each evolved star candidate across multiple wave¬ 
lengths, facilitating the inspection of the local 
stellar and interstellar environment. 

YSOs and dust clumps tend to be concen 


trate d along bright ISM filaments (|Kim et al 


20101 : Sewilo et al. 2013h : conversely, we expect 
there to be very few AGB stars associated with 
interstellar gas or dust structures. Evolved stars 
are generally isolated and inhabit regions of low 
column density. Inspection of the images indi¬ 
cates that many of the HERITAGE-BMC evolved 
star candidates are embedded within bright knots 
of emission, located along large-scale dusty ISM 
structures, or are concentrated around molecular 
clouds; suggesting that these evolved star candi¬ 
dates are instead likely to be YSOs or massive 
stars. We remove sources located within sites of 
massive star formation from our list if their SED 
also resembles that of a YSO. 

Visual inspection of the Spitzer images reveal a 
cluster-like morphology at the location of several 
of our evolved star candidates. This morphology is 
typical of small proto clusters, which appear as a 
single massive object in Herschel. The spatial res¬ 
olution of the HERITAGE images is > 2 pc at the 
distance of the Magellanic Glouds, and blending of 
clustered sources is prevalent. The BMG contains 
a GroupID which specifies group membership and 
is assigned to sources that have multiple matches 
in the BMG. As the majority of the sources in the 
catalog are not sorted into groups we use this cri¬ 
terion as an additional indication of evolutionary 
status and ensure that only isolated sources are 
retained. 

3.3.2. Background galaxies/AGN 

Some nearby galaxies are resolved with Spitzer 
IRAG and can be identified from their morphol¬ 
ogy; for one source in particular (HSOBMHER- 
IGG J85.622736-71.433709) we see a clear spiral 
pattern in the image. Other galaxies (e.g. HSOBMHER- 
IGG J82.644586-68.104609) appear elongated and 
have an edge-on disc morphology, we remove fif¬ 
teen spatially resolved galaxies from the sample. 

As a secondary di agnos tic for background 
galaxies, ISeale et ^ ( 2014 ) identified distant 


galaxies in the BMC via their shape in the SPIRE 
images. Distant galaxies will be isolated from 
interstellar gas and dust structures, and ap¬ 
pear unresolved and thus point-like at the Her¬ 
schel wavelengths. Whereas, other HERITAGE 
source’s, such as YSOs and dust clumps will ap¬ 
pear slightly extended at SPIRE wave l ength s. For 
each source in the BMC, Seale et ^ ( 2014 ) com¬ 
pares the sources 250 /rm full-width half-maximum 
(FWHM) to that of the point spread function 
(PSF). If the values are consistent to within 10%, 
the source is defined to be point-like. Twenty eight 
sources in our sample fulfill this criterion, although 
not all of them are necessarily background galax¬ 
ies. Evolved stars without a detached shell will 
also appear point-like, thus we regard the shape 
in the far-IR as a rather weak constraint for sepa¬ 
rating background galaxies from AGB stars in the 
Magellanic Clouds, and do not remove any object 
on this basis. 

3.4. Luminosity checks 

To further refine our evolved star candidate 
list, we calculate the bolometric luminosity of the 
sources via a trapezoidal integratio n of the SED 


from the optical to the far-IR ('e.g. IWoods et al 


2011 : Jones et al.l 20141: Ruffle et aLT 2015 ). All 
sources within each Magellanic Cloud are consid¬ 
ered to be at approximately the same distance and 
we use a distance modulus of 18.5 to the LMC and 
18.9 to the SMC ( Keller fc Wood 20061) . Where 
possible, we also compute luminosities using the 
IRSF photometry rather than the 2MASS data 
and fold in the photometry from the WISE and 
AKARI surveys. 

Figure [8] shows the absolute bolometric magni¬ 
tudes distribution of the evolved star candidates in 
the LMC and SMC. The luminosites for the ma¬ 
jority of the objects in our sample are consistent 
with evolved stars. Objects more luminous than 
the classical AGB luminosity limit of Mboi = —7.1 
( Wood et al.lll983 ) may be RSGs or massive AGB 
stars undergoing hot bottom burning. At the low- 
luminosity end objects with Mboi > —3.9 will 
be fainter than the tip of the Red Giant Branch 
(RGB) and are unlikely to be an evolved star, as 
stars on the RGB or early AGB do not experience 
significant dust formation. In general these lim¬ 
its are not an absolute division between the RGB, 
AGB and RSG classes, instead they provide useful 
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guidance on the likely class. 

3.5. Literature classification 



Bolometric magnitude 

Fig. 8.— The luminosity distribution of the 
evolved star candidates in the LMC (solid grey) 
and SMC (black-striped). The dashed lines mark 
the expected luminosity range of AGB stars. 


To aid in the final classification, we performed 
a SIMBAeH search on all of the sources to retrieve 
any other information pertinent to the classifica¬ 
tion. We found that twelve of our sources are rela¬ 
tively unknown in the literature, 26 are identified 
as evolved stars or PNe with a high degree of con¬ 
fidence and 5 were classified as a YSO (see Sec¬ 
tion a. We removed known early type stars from 
our list if they were identified by means other than 
colour-magnitude cuts. 


3.6. YSO-fitter screen 


Fina lly we fit the evolved star c andidates us¬ 
ing the Robitaille et al. ( 20061 2007ll grid of YSO 
models, and visually examine the SED and model 
fits. We define well-fitted models as those whose 
per data point is less than ten. 

The YSO-fitter has been shown to be a good 
method of separating evolved stars from YSO can¬ 
didates, since the YSO models do not span the 
high luminosities and r elatively cool photosp heres 
typical of evolyed stars (jWhitnev et al.ll2nn8[l . For 
other populations such as PNe the YSO fitter is 
less effectiye as a discriminant, as YSOs and PNe 
haye similar SEDs and colours. For these sources a 
good fit does not ineyitably imply that the source 
is a YSO. For the evolyed star candidates 3 of 
the 38 sources were well fit by the YSO models; 
these sources are more likely to be true YSOs than 
a post-main sequence star, and thus they are re- 
moyed from the sample. 


4. Results 

4.1. Classification Summary 


From the 529 evolyed stars candidates proyi- 
sionally identified in the HERITAGE BMCs, we 
find strong eyidence for 35 post main sequence 
stars. Eighteen of these sources were confirmed 
spectroscopically while the remaining seyenteen 
sources were selected using the stringent criteria 
described in Section [3] In total 206 objects were 
rejected from our initial BMC sample as their 
SEDs and location suggests these are YSOs; a 


^ http://simbad.u-strasbg.fr/simbad/ 
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small number (< 17) of these sources may poten¬ 
tially be evolved. 

Tables [T] and [2] list our final catalog of evolved 
stars, PNe and SNR with a Herschel detection in 
the LMC and SMC, respectively; their spatial dis¬ 
tribution on the sky is displayed in Figure [9l over¬ 
laid upon the HERITAGE 250 /rm images of the 
LMC and SMC. 

Figures 131 and [7] present the SEDs spanning op¬ 
tical to far-IR wavelengths of each object in the 
sample, and the photometry used to construct the 
SEDs is summarized in Table [3] Finally, Table S] 
provides a tally of the various classes of post main 
sequence stars in the HERITAGE BMC. The fol¬ 
lowing subsections discuss the individual objects 
in our final sample grouped according to their ini¬ 
tial mass and evolutionary status. 


4.2. Comparisons to colour-magnitude di¬ 
agrams 


IR colour-magnitude diagrams have been used 
to classify large numbers of objects detected by 


photometric surveys in the MCs le.g. iBlum et al 


20061 : iMatsuura et al. 2009t iBover et^ 2011 b 


We compare the evolved sources with a Herschel 
detection to mid-IR colour-magnitude diagrams 
of the total evolved stellar population of the LMC 
and SMC. Figure ITOl shows the objects identified 
as candidate evolved stars detected with Herschel 
plotted as blue diamonds in the [8.0] versus [3.6]- 
[8.0] CMD; Figureinishows them in the [24] versus 
[8.0]-[24] diagram. The stars in our sample tend 
to be offset/isolated from the majority of evolved 
star candidates identified from the SAGE sample; 
they are much redder in colour and have bright 
[8] and [24] magnitudes. Due to their very dusty 
nature they typically lie outside CMD cuts for 
evolved stars or reside in areas of high confusion. 

In the LMC there are four exceptions to this. 
These sources are identified as HSOBMHERICC 


J78.424908-67.403028 (BSDL 923), J78.608577- 
68.60194 (2MASS J05142617-6836085), J82.673312- 
69.193322 (HV 2605) and J85.340519-69.529929 
(HV 5999) and are classified as either 0-AGB/PAGB 
stars or blue supergiants. These stars are oxygen- 
rich and have SEDs that show a hot stellar compo¬ 
nent and a second peak at longer wavelengths due 
to a detached dust shell. Thus their [3.6]-[8.0] 
colour will resemble an evolved star with little 



.. 

40"* 20"* 2'‘00’" 40™ 20™ l^OO™ 40™ 0'‘20™ 


RA (J2000) 


Fig. 9.— The BMC post-main sequence sources 
distributed on the sky, overlaid upon the HER¬ 
ITAGE SPIRE 250 /rm images of the LMG (top) 
and SMG (bottom). 
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Table 1: Details of the final sample of evolved star candidates in the LMC. For each source, the BMC name, 
BMC coordinates, bolometric magnitudes (Mboi), total luminosity, SAGE-Spec identification (SSID) and 
final classification is given (see Section |T] for details). 


BMC Source Name 

Alias 

RA 

(J2000) 

Decl. 

(J2000) 

Alboi 

Luminosity 

(Lq) 

SSID 

Classification 

J73.375205-69.296953 

Brey 3a 

04:53:30.0 

-69:17:49 

-5.9 

17784 

4060 

WR 

J73.793532-68.341596 

WOH G064 

04:55:10.4 

-68:20:30 

-9.3 

403277 

4091 

RSG 

J74.195411-69.840206 

R66 

04:56:46.9 

-69:50:25 

-7.1 

52908 

4113 

B[e] 

J74.728949-68.844596 

RP 1805 

04:58:54.9 

-68:50:41 

-1.5 

331 

24 

PN 

J75.5321-71.337365 

R71 

05:02:07.0 

-71:20:13 

-7.1 

55722 

4162 

LBV 

J75.879329-69.967009 


05:03:31.0 

-69:58:01 

-4.8 

6668 


Possible evol. 

J75.924371-71.11932 


05:03:41.8 

-71:07:10 

-5.2 

9874 


Possible evol. 

J76.216758-68.652343 

SMP LMC 21 

05:04:52.0 

-68:39:08 

-3.3 

1715 

4196 

OPN 

J76.989649-68.863313 

SMP LMC 28 

05:07:57.5 

-68:51:48 

-4.1 

3518 

4209 

OPN 

J77.014627-68.671229 

SMP LMC 29 

05:08:03.5 

-68:40:16 

-3.8 

2516 


PN 

J78.424908-67.403028 

BSDL 923 

05:13:41.9 

-67:24:11 

-3.5 

1913 

69 

BSG 

J78.608577-68.60194 


05:14:26.0 

-68:36:07 

-4.6 

5374 


0-AGB/PAGB 

J78.773447-71.652753 


05:15:05.6 

-71:39:10 

-3.4 

1776 


Possible evol. 

J78.854807-66.317713 


05:15:25.0 

-66:19:04 

-3.4 

1825 


Possible evol. 

J79.075912-71.899396 

IRAS 05170-7156 

05:16:18.2 

-71:53:58 

-5.4 

11702 

78 

UNK 

J79.617143-68.067639 

IRAS 05185-6806 

05:18:28.1 

-68:04:04 

-4.9 

7201 

4307 

G-PAGB 

J81.264647-69.642784 

N 132D 

05:25:02.0 

-69:38:39 




SNR 

J81.920736-69.135316 

IRAS 05280-6910 

05:27:46.2 

-69:08:07 

-7.8 

107347 

4454 

RSG 

J82.673312-69.193322 

[BE74] 328 

05:30:41.5 

-69:11:36 

-4.7 

6224 


B[e] 

J82.683128-71.716888 

IRAS 05315-7145 

05:30:43.9 

-71:43:01 

-5.0 

7688 

125 

C-AGB 

J82.730596-68.59895 


05:30:55.3 

-68:35:56 

-4.2 

3938 


Possible evol. 

J83.558492-69.789067 

IRAS 05346-6949 

05:34:14.0 

-69:47:21 

-7.7 

91262 

4555 

RSG 

J83.866283-69.269995 

SN 1987A 

05:35:28.0 

-69:16:11 




SNR 

J84.10768-69.381998 

R126 

05:36:25.8 

-69:22:55 

-7.3 

64377 

4601 

B[e] 

J84.173847-66.808723 


05:36:42.0 

-66:48:31 

-3.8 

2529 


Possible evol. 

J85.156861-70.153736 


05:40:37.6 

-70:09:13 

-3.8 

2543 


Possible evol. 

J85.340519-69.529929 

HN 5999 

05:41:21.7 

-69:31:48 

-7.1 

56957 


0-AGB 

J87.249164-70.556065 

IRAS 05495-7034 

05:48:59.7 

-70:33:22 

-5.3 

10049 

190 

C-AGB 


Note: The HERITAGE-BMC designation prefix is HSOBMHERICC. 


Table 2: Details of the final sample of evolved star candidates in the SMC. For each source, the BMC name, 
alias, BMC coordinates, bolometric magnitudes (Mboi), total luminosity, SMC IRS identification (SMC IRS) 
and final classification is given. 


BMC Source Name 

Alias 

RA 

(J2000) 

Decl. 

(J2000) 

Alboi 

Luminosity 

(Lq) 

SMC 

IRS 

Classification 

J12.035323-73.248808 

LIN 107 

00:48:08.5 

-73:14:56 

-9.0 

324932 


PN 

J12.150026-72.967348 

LIN 115 

00:48:36.0 

-72:58:02 

-6.2 

24197 

32 

CPN 

J12.742224-73.86164 


00:50:58.1 

-73:51:42 

-4.6 

5427 


Possible evol. 

J12.903066-73.338318 

Jacoby SMC 17 

00:51:36.7 

-73:20:18 

-3.3 

1624 


B[e] 

J13.050807-73.147819 

BMB-B 75 

00:52:12.2 

-73:08:52 

-5.9 

17561 

177 

0-AGB 

J13.579993-72.485754 


00:54:19.0 

-72:29:09 

-5.2 

9946 

260 

Symbiotic star 

J20.652614-73.171815 


01:22:36.6 

-73:10:19 

-3.7 

2404 


Possible evol. 


Note: The HERITAGE-BMC designation prefix is HSOBMHERICC. 
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Table 3: Summary of the columns present in the table of photometry which is available on-line only. 


Column 

Number 

Column Name 

Description 

Null 

Value 

1 

SourceName 

HERITAGE-BMC Designation 


2 

RA 

Right Ascension (J2000) 


3 

Dec 

Declination (J2000) 


4 

Galaxy 

LMC or SMC 


5-12 

fU, dU, fB,dB, fV, dV, fl, dl 

MCPS U, B, V, I £iuxes (£) and errors (d) [mjy] 

-999.9 

13-18 

fj, dj, fH, dH, fK, dK 

2MASS J, H, K £luxes and errors [mJy] 

-999.9 

19-26 

£3.6, d3.6, £4.5, d4.5, £5.8, d5.8, £8.0, d8.0 

IRAC fluxes and errors [mJy] 

-999.9 

27-30 

£24, d24, £70, d70 

MIPS fluxes and errors [mJy] 

-999.9 

31-40 

£100, dlOO, £160, dl60, £250, d250, £350, d350, £500, d500 

PACS & SPIRE fluxes and errors [mjy] 

-999.9 

41-45 

fiaglOO, flagl60, fiag250, flag350, flag500 

HERITAGE-BMC quality flags 

-999.9 

46-53 

£W1, dWl, W2, dW2, £W3, dW3, £W4, dW4 

WISE flux and errors [mJy] 

-999.9 

54-63 

£N3, dN3, £S7, dS7, £S11, d£Sll, £L15, d£L15, £L22, d£L22 

AKARI flux and errors [mJy] 

-999.9 

64 

SSID/ SMC IRS 

Spectral Identification number of the target from 
Woods et ah ('2011'!: Ruffle et ah ('2015'! and 
Woods et ah (in prep). 

-999.9 


circumstellar excess, while at [8.0]-[24] the contri¬ 
bution from circumstellar dust becomes significant 
and their colours resemble the other stars in our 
sample with a strong dust excess. 


4.3. Low and intermediate mass stars 


Table 4: Classification groups and tally for the 


LMC and SMC. 


LMC 

Type 

Count 

SMC 

Type 

Gount 

B[e] 

3 

B[e] 

1 

BSG 

1 

0-AGB 

1 

C-AGB 

2 

PN 

2 

C-PAGB 

1 

Possible evol. 

2 

LBV 

1 

Symbiotic star 

1 

0-AGB 

1 



0-AGB/PAGB 

1 



PN 

4 



Possible evol. 

7 



RSG 

3 



SNR 

2 



UNK 

1 



WR 

1 



Total 

28 

Total 

7 


Our sample contains thirteen confirmed low to 
intermediate mass stars and PNe, which make up 
32 per cent of the total sample. If we include the 
possible evolved stars (Section 14.6L which have 
unknown mass, this value could rise to 60 per cent 
of the probable evolved sources detected in the 
HERITAGE-BMC. Due to their intrinsic bright¬ 
ness, stars with masses below SMq were not ex¬ 
pected to be recovered by the HERITAGE survey, 
as only dust ejected by a few of the most massive 
stars were predicted to have luminosities greater 
than the Herschel sensitivity limits. 


The intermediate mass stars we detect include 
two of the most extreme AGB stars in the LMC 
with exceptionally high mass-loss rates, whose car¬ 
bon rich c hemistry has been co nfirmed spectro¬ 
scopically ( Gruendl et al.l l2008ll . The SEDs of 
the two carbon stars (C-AGB) and the symbiotic 
(AGB-white dwarf binary) star do not show any 
evidence for excess emission in the far-IR, how¬ 
ever. All the oxygen-rich AGB (0-AGB) stars in 
our sample show a stellar component, a compo¬ 
nent due to warm dust and in some instances a 
far-IR excess may also be present. This excess may 
be associated with the dust-driven wind, a stable 
circumstellar disc or dust swept-up from the sur¬ 
rounding ISM. The mass-loss rate stated for the 
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Fig. 10.— The [8.0] vs. [3.6]-[8.0] CMD for the 
LMC (top) and SMC (bottom) showing the dis¬ 
tribution of the post-main sequence stars in our 
final sample (blue diamonds) plotted over a Hess 


diagra m of the evolved star sample of IB over et al 
(1201 ih for the LMC and SMC. 


Fig. 11.— The [24] vs. [8.0]-[24] CMD for the 
LMC (top) and SMC (bottom) showing the dis¬ 
tribution of the post-main sequence stars in our 
hnal sample (blue diamonds) plotted over a Hess 


diagra m of the evolved star sample of IB over et al 
(120111) for the LMC and SMC. 


15 










0-AGB stars is likely an underestimate as cur¬ 
rent model grids do not account for complex dust 
geometries due to a detached shell. Detailed 3D 
radiative-transfer modeling of the sources with a 
far-IR excess is required to determine reliable dust 
production rates and to differentiate between ISM 
dust and dust produced in the stellar wind (Jones 
et ah, in preparation, Paper II). 

Our final sample also includes several highly 
evolved post-superwind objects such as planetary 
nebulae. The five planetary nebulae in our sample 
have been identified from a combination of Ha im¬ 
ages, optical ionized gas emission lines and mid-IR 
spectra. Due to the similarities between the SEDs 
of PNe and YSOs and the strict criteria we ap¬ 
plied when selecting the evolved star sample it is 
likely that we have excluded some PNe from our 
final list. Thus the number of PNe with a Herschel 
counterpart included in our final catalog should be 
treated as a lower limit. 


4-3.1. AGB stars 

HSOBMHERICC 378.608577-68.60194 (2MASS 
J05142617-6 836085): This lon g period variable 
(P = 382 d; Fraser et al. l2008ll is likely oxygen- 
rich. It was found to ha ve a mass-loss rate (MLR) 
of 2.8 X 10"^ Mgyr-i (|Riebel et al.l 1201211 . The 
SED is double peaked which may indicate that the 
source is a post-AGB star or that it has an unre¬ 
solved detached shell. 

HSOBMHERICC J82.683128-71.716888 (IRAS 
05315-7145) is an extreme carbon-rich AGB star, 
heavily obscured by dust. It has a v ery high mass- 
loss rate of 1.7 x 10“"^ M 0 yr“^ ( Gruendl et al 


I 2 OO 8 II : this is the second highest of all the extreme 
AGB (x-AGB) stars in the EMC. SiG absorption 
features are present in the IRS spectra. 
HSOBMHERICC J85.340519-69.529929 (HV5999): 
This source has a luminosity of Mboi = —7.1 
mag and a spectral type of M3 suggesting the 
source is likely a massive 0-rich AGB star un¬ 
dergoing hot bottom burning or is a red super- 
giant. It is also includ ed in the RSG catalog of 
liviassev fc Olsen ( 2003ll . who give Mboi = —8.05 
mag. Fits with GRAMS indicate that this source 
is oxygen -rich with a mass-l oss rate of 1.4 x 10“® 
MQyr“^ ( Rlebel et al. 2012 1. The SED shows ey- 
idence of a far-IR excess which (if associated with 
the star) may indicate the presence of large (^m 


- mm size) grains and/or cold dust in a detached 
shell. 


HSOBMHERICC J87.249164-70.556065 (IRAS 
05495 - 7034 ) is an extremely red ca rbon star first 
identified by iGruendl et al. ( 2008h . Its spectra 
show the presence of amorphous carbon, SiG ab¬ 
sorption features, and a strong MgS feature at 30 
p,m. This ‘extremely red object’ (ERO) has the 
highest mass-loss rate (2.3 x 10~' ^ M(7iyr~ ^) of 


all the x-AGB stars in the LMC (iGruendl et al 

I 2 OO 8 II . 

HSOBMHERICC J13.050807- 73 .14 7819 (BMB-B 
75) is a Mir a-type variable with a long period of 
P = 1453 d (iGioni &: Habin g 2003ll . It is an M 6 - 
type star (IBlanco et al.lll980ll . The star is oxygen- 
rich with a mass-l oss rate of 1.5 x 10“® M 0 yr“^ 
(IJones et al.ll2012h . Recently a redshifted galaxy 
(z~0.16) has been discovered in the same line of 
sight as BMB-B 75 (Kraemer et al. 2015, in prep.) 
The far-IR emission detected by Herschel is prob¬ 
ably due to this spatially coincident galaxy. 
HSOBMHERICC 313.579993-72-485754 (SSTIS- 
ACEMA 3005419.21-722909.7) is a D-type sym¬ 
biotic star, consisting of an AGB star and white 
dwarf (|01iveira et al.N2013r) . Of the seven known 
symbiotic stars in the SMG it is the only one with 
warm dust; this dust is silicate rich and unpro¬ 
cessed as it shows no PAH emission. 


4 . 3 . 2 . Post-ACB stars 


Our sample only inc ludes one confirme d post- 
AGB star. However ISeale et al.l (j2014l l notes 
that the B MG contains four confir med post-AGB 
stars from Ivan Aarle et al.l (l201ll) with low res¬ 
olution optical spectra. Of these four sources, 
both HSOBMHERIGG J73.375205-69.296953 and 


J78.854807-66.317713 are included in our sam¬ 
ple as a WR star and a probable evolved star, 
respectively (see Sections 14.41 and 14.61) . How¬ 
ever, HSOBMHERICC J80.919272-68.090965 and 
J83.220747-69.987453 were eliminated from our 
sample either because the SED resembles that of 
a YSO or because the SED at longer wavelengths 
appears disjointed and the Herschel data were at¬ 
tributed to a mismatch. The SEDs of these two 
sources are shown in Figure [2j 


HSOBMHERICC 379.617143-68.067639 (IRAS 
05185-6806) w as first identified as a likely post- 
AGB object by yan Loon et al. 1 1997ll . Its mid- 
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IR Spitzer-IRS spectrum shows the presence of 
the unidentified 21-/im featur e found only 


carbon-rich post-AGB stars (fVolk et al~ 2011 


Matsuura et al. 2Ql4 ISloan et ah 2014ll . confirm- 
ing that mass loss during the AGB phase has 
terminated. 

4.3.3. PNe 

HSOBMHERICC 376.216758-68.652343 (SMP 
LMC 21) \s a. very high-excitation, qua drupolar 
PN of Type 1 ( Stanghellini et al1ll999ll . It has 
an extremely strong IR-excess (Tdustcont = 130 
K) and the Spitzer IR spectrum shows promi¬ 
nent crystalline silicate emission due to forsterite 
and enstatite. No other sol id-state features are 
seen ( Stanghellini et al. IH). The spectrum also 
shows both low and very high excitation nebular 
emission lines (i.e., [S ill] and [Ne vi], respec¬ 
tively). It is one of the reddest PNe in the Mag¬ 
ellanic Clouds. The ionized mass in the n ebula is 
determined to be between 0.2 - 0.6 Mq ( Barlow! 
19871 : iBoffi fc Stanghellinilll994l . and the progen¬ 
itor is thought to be a relatively massive AGB 
star. 

HSOBMHERICC J76.98964 9- 68.863313 (SMP 
LMC 28) \s a high-excitation isolated, oxygen- 
rich PN (|Reid__fc_P arkei 20061 : Hora et ^ I 2 OO 8 I : 


Bernard-Salas et al 


20m 


It is has a high ni- 
troge n and neon abundance compared to other 


PNe ( Stanghellini et al.lbOOSt Bernard-Salas et al 


20081) . This PN is resolved in HST images 
and has three co-aligned round shapes, corre- 
sponding to the c entral star and two faint arms 


(Shaw et al. 200 


bv iVillaver et al 


]) . The central star was detected 


( 2 OO 3 II . They estimate its lumi¬ 


nosity to be L = 1.4 X IO^Lq, with Teg < 17000 
K, which corresponds to a progenitor mass of 
^ 3M0. The ionized mass in the nebula is 


< O. 2 M 0 (IWood et al .111 98711 . 

HSOBMHERICC J77.014627-68.671229 (LHA 
120-N 102/ SMP LMC 29) i s a bright PN no ted to 
have strong emission lines ( Reid fc Parked [ 2 OO 6 II . 
This PN is resolved and has an irregular shape 
with a bipolar core in HST images. It is thought 
to be one of the smallest bipolar PNe recog- 
nise d to date, with ra dius of ~0.27 arcsec in [O 
III] ( Shaw et al. 20061) . No central star is evi¬ 


dent within the brig ht nebulosity of the image 
(Villaver et al. 20031). It has emission lines in 


extre me photoionization tem peratures (T 
kK) ( Herald fc Bianchil 2007 ). 


200 


HSOBMHERICC J12.150026-72.967348 (Lin 115/ 
SMP SMC 11) is a well known carbon-rich PN 
which was first discovered by Henize ( 19561) . Both 
emission line observations and PAH features in 
its IRS spectrum c onfirm its carbon-rich 


na¬ 


ture (lStarigh£]lini_et_alJ 12003 |BemMd^Sa]as_et_al 
2 OO 9 I) . High-resolution HST images show that 


it has a complex bipolar morp hology, which ex¬ 


tends out to 2" from the center (|Stanghellini et al 
200,31 ). It has strong mid-lR emission; the dust 


continuum sharply rises after 12 /rm and its SED 
peaks around 35 rim, indic ative of a far-lR-excess 
( Bernard-Salas et al.ll2009l) . 

HSOBMHERICC J12.035323-73.248808 (LIN 
107) was identified as a planetary nebula b y 


Lindsay ( 1961 ): Mevssonnier fc Azzonardil ( 199.3f ). 

Other authors list this object as a weak emission- _ 

line object or ‘low-excitat i on blob’ (iMevnadier &: Hevdari-Malaveri 
2007 : Sheets et al. 2013t iTestor et al. 20141) sug- 
gesting that it is an embedded star cluster or 
compact Hii region. 

4.4. Massive stars 

We have compared the HERITAGE-BMG to 
the Spitzer ca talogs of massive sta r s in t he Magel¬ 
lanic Glouds ( Bonanos et al. 20091 20101 ) and find 
nine matches in the LMC and one match in the 
SMC. This is approximately twice the number an¬ 
ticipated prior to the HERITAGE observations. 

Notably far-lR emission is observed for four of 
the fifteen confirmed B[e] stars in the Magellanic 
Clouds, while only one of the 106 late-type WR 
stars has a far-lR detection. 

Three heavily dust enshrouded extreme RSGs 
are also detected. These stars all optically thick 
in the mid-lR where amorphous silicate absorp¬ 
tion features are seen at 10 and (in some in¬ 
stances) 20 /im indicative of extremely high mass- 
l oss rates, which may be as large as 10“^ M 0 yr“^ 


(jDiikstra et al.N2003l) . The mass loss from these 


three sources w as missed in dust bu dget estimates 
for the LMC bv iRiebel et al. ( 20121) . 


high states of ionization (e.g., N v) which have 


4.4-1. BSC stars 

HSOBMHERICC J78.424908-67.403028 (BSDL 
923) is a blue supergiant (BSC); a massive star, 
newly evolved from the main sequence. Its loca- 
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tion in the young cluste r LMC-N30 (|Gouliermis et al 
2003; Bica et al. lOOOl l is consistent with that of 
a massive star. The SED is distinctly double- 
peaked, falling to a minimum at '^8 /im before 
rising steeply long-ward of 15 ^m; this red peak 
is due to a circumstellar disc or torus. Weak sili- 
cate features are se en in its Spitzer-IRS spectrum 
( Woods et al. 2011 1. 


4-4-2. RSG stars 

HSOBMHERICC J73.793532-68.341596 (WOH 
G064 ^ is a very luminous red supe rgiant with spec¬ 
tral type M7.5 (jElias et al.1119861 1 and an effective 
temperature of Teff 3200 K. It is surrounded by 
a thick dust shell and is one of six evolved stars in 
the Magellanic Clouds to have a 10 g m amorphous 


silica t e feature in self-absorption (lRoch eet_al 


1993 : Warns et ahl 19991 : Buchanan et al.l 120061 : 


Kemper et al.l 1201(1 Jones et al. (in prep)). Its 
mass-loss rate is estimated to be in th e region 
of 10“^ M 0 yr“^ ( van Loon et al. 2005all with a 
wind speed of Uwind ~ 25kms~^, dete rmined from 
OH, SiO, and H:;)0 maser emissi o n (fu^od et al 


Ivan Loon et al.l Il996l 119981 : iMarshall et al 
2004ll. The circu mstellar dust shell was resolved 


by lOhnaka et al.l (|2008r) , who found that the en¬ 
velope had a bipolar geometry; they derived a 
circumstellar envelope mass of 3-9 Mq. 

We estimate the bolometric luminosity to be 
4x 10® Lq (Mboi = —9.3), however if the geometry 
of the dust shell is non spherical then the bolomet¬ 
ric luminosity is an overestimate. Accounting for 


the no n spherical nature of the shell lOhnaka et al 


( 2008l l determine the luminosity to be 2.8 x 10® Lq 
(Mboi = — 8 . 8 ) leading to a down-ward revision of 
the star’s initial mass to Mzams ~ 25M0. Its loca¬ 
tion beyond the Hayashi limit in the H-R diagram 
indicates th at the star is no longe r in hydrostatic 
equilibrium (jLevesoue et al.l 120071 1 and it may be 
experienci ng an intense un s table p hase of heavy 


e pr 

mass loss. Ivan Loon et al.l (l2010bll suggests that 


the mass-loss rate of WOH G064 has increased 
dramatically in the past few thousand years (i.e. it 
has experienced the onset of the superwind phase), 
as they find no evidence for a submm excess. 
HSOBMHERIGG J81.920736-69.135316 (IRAS 
05280-6910) is a highly obscured OH/IR red su- 
pergiant with an initial mass Mzams 2OM0 
(|Wood et 'n] I 1992 I) . Its stellar counterpart was 
identified using high-resolution near- and mid-lR 


imagery by Ivan Loon et al.l (l2005b[l . OH maser 
emission has been used to determi ne its wind out¬ 
flow velocity which is ss 20 kms~^ ( Marshall et al 


2004); the star also exhibits H 2 O maser emission 


(Ivan Loon et al.l 1200 111 . It has an extremely high 
mass-loss rate of M cs 2 — 8 x 10“® M© yr“^ 
( Bover et al. 2010l ) and is potentially the most 
dust-enshrouded supergiant in the Magellanic 
Clouds. In the IRS spectra both the 10 and 20 
pm amorphous silicate features are seen in absorp¬ 
tion and a declining, fe atureless continuum is seen 
at far-IR wavelengths ( Kemper et al. 2010l) . The 
dust shell may have a flattened geometry which 
(unlike WOH G064) is viewed edge-on. The total 
dust mass is « 0.3 M 0 . There is no evidence of ex¬ 
cess emission at A > 100 pm, implying that there 
is no significant contribution from dust colder than 
'^lOO K or fro m large (um - mm size) grains in a 
detac hed shell ( van Loon et ahlEoiOb : Bover et al 
20101 ) . 


HSOBMHERIGG J83.558492-69.789067 (IRAS 
05346-6949): This luminous oxygen-rich red su - 
pergiant is highly obscured ( Elias et ^ Il986[l : 
due to its very red IRAC [3.6]-[4.5] ^ 2.5 
colour it has been classified as an ERO by 
iGruendl &: Chu (12009) and an extreme-AGB star 


by Srinivasan et al. ( 20091) . We classify it as a 


RSG because of its high luminosity an d the silicate 


features present in its IRS spectrum () Jones et al 


20141 ). Unlike many sources in the ERO class 


this source is not carbon-rich, instead it has an 
oxygen-rich dust composition. Like WOH G064, 
the Spitzer-IRS spectrum of HSOBMHERIGG 
J83.558492-69.78906 7 exhibits a 10 p m silicate 
feature in absorption ( Jones et al.ll2012 . Woods et 
al. (in prep), Jones et al. (in prep)) which means 
that the dust envelope is optically thick, and the 
star is experiencing heavy mass loss. Despite be¬ 
ing bright in the m id-IR, any maser emi ssion has 
not been detected (jyan Loon et al.l 1200 ij) and the 
source is not bright in the optical. Radio con¬ 
tinuum emission at 1420 MHz may be associated 
with the source. 

4.4-3. B[e] stars 

HSOBMHERIGG J74.195411-69.840206 (R66) is 
a well-known B[e] supergiant of spectral type B 8 
la, with luminosity L = 3 x 10®, Tes ~ 12000 
K and an in i tial mass of Mzams ^ 3 OM 0 
( Lamers et al.l Il998t Stahl et al. 1983 ). The 
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star is both coo l and below t he H umphreys- 
Davidson limit ( Zickgraf et al. IQSfih . leading 


van Loon et al. ( 2010bll to suggest that it is a blue- 


loop star on its way back to becoming an RSG. 
The bpjtzer-IRS spectrum shows PAHs, amor¬ 
phous silicates and strong crystalline silicate fea¬ 
tures, indicative of dust grain proce ssing in a long- 
lived pole-on circumstellar disc (jKastner et al 
2 OOIIL 2010 ! : Aret et al. 2012tl . This disc is mas¬ 
sive and contains ss 0.001 Mq of dust. 
HSOBMHERICC J82.673312-69.193322 ([BElj] 


328) is a luminous (Mboi = —10.15; iMassev 


20021). emission line star of spectral type B1.5IIe 
( Reid fc ParkeiJ20l^ . It is listed as a poorly stud¬ 
ied eruptive variable of early spectra l type in the 


Gene ral catalog of Variable Stars (iSamus et al 
20041) . 


HSOBMHERICC J 84 .10768-69.381998 (R126) is 
an extremely luminous B[e] supergiant in the 
LMG. Its bolometric luminosity is L = 1.2 x 
1 O®L0, with Teff ^ 22500 K and is estimated 
to have an initial mass of Mzams 75Mr:^. Its 
prese nt day mass is M ^ 40Mq ( Zickgraf et al. 
19851) . Estimates for the current mass-loss rate 
range between M = 10“® — 4.6 x 10“® Mq 
yr“^, depending on the assumed terminal veloc- 
i ty of the stellar wind Vnn ~ 650 — 1800 km s“^ 


( Zickgraf et al.l [l985 : Biorkman et al. 19981) . R 


126 is thought to have a relatively massive ‘puffed 
up’ circumstellar disc which is estimated to con¬ 
tain Mdust 
20061 ) . 


0.003 Mq of dust (jKastner et al 


Detailed analy s is of R 126’ s disc can be 

(l2007l) . 


found in iPorteil (120031 ) ; iKrans et al 


HSOBMHERICC J12.903066-73.338318 (Jacoby 
SMC 17) is included in the list of probable 
B[el supergiants in the SMG ( Mennickent et aP 


20021 ) . as its light curve resembles Gala ctic B[e 


stars. Mevssonnier fc Azzopardil ( 1993h identify 


the source as a very low excitation object with a 
blue continuum. It was also identif i ed as an emis¬ 
sion line star bv lMurp hy fc Besselll (2000 ). It may 
potentially be a PN ( Boroson fc Lieberti I989I) ; its 
total luminosity is on the low-side for a B[e] star 
but is consistent with a PNe. 

4.4-4- LBV stars 

HSOBMHERICC J75.5321-71.337365 (R71) is 
one of six confirmed dust y luminous blue variables 
(LBVs) in the LMG ( Humphreys fc DavidsonI 


19941) . Its initial mass is esti mated to be Mzams 
40Mq ( Lennon et al. 1993 ). LBVs can develop 
instabilities and show sudden bursts of mass 
loss; this mass-loss history can be probed via 
dusty ring nebulae often associated with the star 
( HutsemekersI 1994h . R71 recently experienced 

an eruption, where there was a ^2 mag increas e 


in the V-band brightness ( Szczvgiel et al. 20101) . 


This increase in emission is not seen in our mid-IR 
data which show a smooth transition between the 
Spitzer and Herschel photometry, even-though the 
Herschel data were obtained post-eruption, while 
the Spitzer observations were obtained prior to 
the outburst. The difference in flux between the 
WISE bands at 3.4 and 4.6 pm and the IRAG 
bands may be due to the increased emission from 
the photosphere during the recent outburst. 

The dust shell of R 71 contains PAHs, and both 
amorphous and crystalline silicates indicative of 
non-equilibrium chemistry. It estimated to con¬ 
tain 0.02 Mq of dust, ejected at a rate of 7 x 10“® 
Mq yr“^ over the past 3000 years. A second dust 
shell from an earlier exp ulsion event is esti mated 
to contain 0.3 Mq of dust ( Voors et al ]|l999l) . This 
se cond dust shell has bee n ca lled in to question 


bvivan Loon et al.l (l2010bD and iGuha Nivogi et al 


(|2014D who find that R 71 lacks dust colder than 
Td < 100 K. 

Using data from t he HERITAGE S cience 
Demonstration Program Bover et al. ( 2010l) mea¬ 
sured a far-IR dust component at A > 250 pm. In 
later HERITAGE photometric data-sets, which 
better accounts for background emission, the 
SPIRE 350 and 500 pm is conf irmed to be diffuse. 
Using this refined photometry, iGnha Nivogi et ^ 


(I 2 OI 4 I) have fit radiative transfer models to R71 
and derive a total dust mass of 0.01 Mq and a 
time-averaged dust mass-loss rate of 2.5 x 10“® 
Mq yr-1. 

4.4-5. WR stars 

HSOBMHERICC J73.375205-69.296953 (Brey 
3a/ BAT99-4 ) is & WC9-type Wolf-Rayet (WR) 
star with a prog enitor mass believed to be in ex¬ 
cess of 25 Mq ( Brevsacher et al.l 19991) . Late- 
type WR stars can lose mass a t rates of up to 
10“'^ MQyr“^ ( Growthei 2007). Characterised 
by strong He, N and C emission lines and red 
SEDs due to warm circumstellar dust, WR stars 
are relatively easy to identify; there are thought 
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to be 134 WR stars in the LMC. Brey 3a is un¬ 
usual as its exact spectral type remains contro¬ 
versial, leading some authors to dispute its WR 
status and instead suggest that it may be a tran- 
sition object between an Of star and a WR sta r 
dMoffat 1991 : He vdari-Malaveri fc Melnick 1992 1. 
Egan et al. ( 200 ill note the presence of an M3-type 
red giant 5.2" away from Brey 3a, which may con¬ 


tribute some flux at <24 /rm but Ivan Loon et al 


(|2010bll find that the far-IR emission seen in the 
MIPS-SED spectra of Brey 3a is typical of a WR- 
type star rather than an M-type giant with warm 
circumstellar dust. The IRS spectrum, which rises 
steeply at ^15 gm, also resem bles that of the WR- 
type star ([Woods et al.ll201l[l . 


4.5. Supernova remnants 

Our Herschel sample includes two well known 
supernova remnants: 

A detailed analysis of SN 1987A (HSOBMHER- 
ICC J83.866283-69 269995) in the far-IR/sub-mm 
was conducted by iMatsuiira et al.l (l201ll l2015ll : 


Indebetouw et al. ( 2014ll . who found that the 


SN ejecta contained a large reservoir of freshly- 
synthesised cold dust grains with a temperature 
of ~20 K and a dust mass of 0.4-0.8 Mq. This 
dust mass (measured 24 years after the explo¬ 
sion) is significantly higher than the IO^'^Mq of 
dust reported in the first couple of y ears follow¬ 
ing the explosion ( Wooden et al. 1993ll : providing 
the first clear evidence that supernovae can pro¬ 
duce significant amounts of dust and that they 
can be an important source of dust in the ISM 
of galaxies. The dust produced in the ejecta is 
thought to be composed of silicates and amor¬ 
phous carbon. As the core-collapse SN is rapidly 
evolving and the expanding ejecta have yet to in¬ 
teract with the warm (~170 K) circumstellar dust 
ring (ejected from its re d supergiant progenitor; 
Barkat fc Wheeler 1988h and experience reverse 
shocks which may potentially destroy or alter 
the gr ains (e.g..ljones et al.lll994 iLakicevic et al 


201,4 iTemim et al.l 20151) . the total quantity of 


dust ejected into the ISM by SNe remains uncer¬ 
tain. 

HSOBMHERICC J81.264647-69.642784 (N132D) 
is the remnant of a core-c ollapse supernova tha t 
exploded '^2700 years ago ( Williams et al 
The supernova remnant was detected as a slightly 
extended point source in two Herschel bands. Un¬ 


like SN 1987A, the emission from N 132D at these 
wavelengths is too faint to accurately constrain 
the dust mass of the ejecta. N132D is one of 
the brightest X-ray remnants in the Magellanic 
Clouds, and its extended emission has been ex- 
tensively studied i n optical and IR wavelengths 
(iMorse et al.l 119951 ): the bright mid-IR emission 
has been attributed to swept-up ISM dust colli- 
sionally heated by the hot X-ray e mitting plasma 
gener ated by SN blast waves ( Williams et al.l 
120061 ) . Spitzer observations also reveal large poly¬ 
cyclic aromatic hydrocarbons (PAHs) that have 
been processed by therm al sputtering in the blast 
wave ( Taope et al.ll2006l) . Dust created in th e SN 
ejecta has been observed bv iRho et al. ( 20091) . via 
IR spectral lines; they estimate the mass of the 
newly formed dust in the ejecta to be >0.008 Mq. 

The far-IR emission from othe r Magellanic 


Cloud supernova remnan ts (see lOtsuka et al 


2 OIOI : Lakicevic et al.l[2015 , and references therein) 


are not included in the HERITAGE point source 
catalogs as the dust emission from these rem¬ 
nants is either too faint or it is extended and large 
enough to be resolved by the Herschel observa¬ 
tions. 

4.6. Possible evolved stars 

Two sources in our sample have low signal- 
to-noise or very unusual/unidentifiable features 
in the correspondi ng gpfteer-IR S spec tra. Due 
to this uncertainty IWoods et al.l (I 2 OIII) classified 
these sources as unknown objects. As the true 
nature of these objects has not been spectroscopi¬ 
cally confirmed, we assign these objects a possible 
evolved star status. 

HSOBMHERICC J74.728949-68.844596 is proba¬ 
bly a planetary nebula (RP 1805); it has a diam¬ 
eter in Ha of 5 arcsec ( Reid fc Parker 2006[) . The 
mid-IR IRS-spectrum of the source shows no clear 
identifying features above the noise and was thus 
reported as a non-detection. The source is faint 
and appears elongated the SAGE-LMC 8-/xm im¬ 
ages. A UV-bright star located within 1 arcmin 
of the source could conceivably be irradiating the 
dust cloud. 

HSOBMHERICC J79.075912-71.899396 (IRAS 
05170-7156) is an unusual object. It was as- 
signed an unknow n (UKN) spectral type by 
Woods et al.l ( 2011 ). The mid-IR spectrum rises 
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toward longer wavelengths and is featureless apart 
from a dip in the continuum between 11.5-16.0 
/rm. Its SED is indicative of a post-main se¬ 
quence object, however, it has also been con- 
sidered to be an AGN candidate (Ide Griio et al.l 


1987 : Kozlowski &: KochanekI 2009ll and a hi gh 


probability Stage 1 YSO ( Whitney et al.ll20(M . 

Nine sources in our sample are relatively un¬ 
studied in the literature, but detailed examina¬ 
tion of the SEDs and fits with YSO models sug¬ 
gest that these sources are likely evolved stars. 
Only one source HSOBMHERICC J84.173847- 
66.808723 was given a less certain classification 
due to its unusual multi-peaked SED and other 
characteristics which meant that an alternate clas¬ 
sification could not conclusively be ruled out. 

4.7. Physical extent of dust shells 

Bow shocks due to wind-lSM interactions and 
detached shells around AGB stars are commonly 
seen in the Milky Way via their far-IR excess emis- 
sion dKerschbaum et aDl2010l: iGroenewegen et al.l 


201 ll : IGox et aU 120121 1. High-mass evolved stars 

are often associated with a ring nebula formed 
from either past stellar eruptions or by swept up 
ISM material. In all instances the maximum phys¬ 
ical size observed between the central star and 
the d etached shell or bow shock is smaller than 
1 pc ( Gox et ^l2012ll . These size scales are much 
smaller than the highest achievable resolution of 
8.6-8.8", 2-3 pc at 100 urn in the HERITAGE 
images. 

Given the physical size of the Herschel beam, 
we cannot use the images to distinguish between 
dust formed in the ejecta or dust swept-up from 
the surrounding ISM. Instead, the spectral energy 
distributions may provide some indication of far- 
infrared excess emission. The circumstellar en¬ 
velopes of the extreme carbon stars, RSGs, LBV 
and WR stars in our sample lack a far-IR excess, 
and thus are unlikely to have large (/rm - mm size) 
grains or harbor any significant reservoirs of cold 
dust. However, we may see excess emission (i.e. 
above what is expected from a circumstellar shell) 
from some of the oxygen-rich evolved stars, plan¬ 
etary nebulae and B[e] supergiants. This emission 
may either be due to dust produced by the evolved 
star or it may arise from swept-up ISM material. 

Dust mass estimates obtained using graybody 


fits to the far-IR emission can in principle be used 
to identify the li kely source of this emission (e.g. 
Seale et al. 20141 1. however this method is unlikely 
to yield reliable dust masses due to the complex 
geometry of the objects in our sample. Instead 
3D radiative transfer modeling is required to de¬ 
rive reliable dust masses and temperatures (Jones 
et al. in prep. Paper H). High-resolution follow¬ 
up observations with the Atacama Large Millime¬ 
ter/submillimeter Array (ALMA) should also un¬ 
cover the true source of this far-IR emission. Fol¬ 
lowup observations will probe the mass-loss rate 
history of the star or the wind-ISM interaction at 
low metallicity. 

4.8. Selection effects and completeness 

The evolved stars in our sample were identified 
from the HERITAGE far-IR BMG, as such they 
are subject to the same biases, selection effects 


and li mitation discussed in detail by iSeale et al.l 
(I 2 OI 4 II . Notably, the HERITAGE SPIRE observa¬ 
tions are more sensitive than the PACS observa¬ 
tions, thus faint sources near the sensitivity limit 
are detected only at 250 and 350 /rm. This will ad¬ 
versely affect the detection of evolved stars which 
are expected to be brighter in the PACS bands 
compared to SPIRE. 

In the HERITAGE images complex structures 
of knots and filaments due to ISM dust dominate 
the background emission; this prevents the detec¬ 
tion of faint sources in regions of high background. 
While AGB stars are unlikely to be directly asso¬ 
ciated with dust clumps or filaments, many may 
lie along the same line of sight. This is even more 
problematic for massive stars which are preferen¬ 
tially found in higher ISM column densities. High 
background {B 250 > 10 MJy sr“^) regions cover 
16% of the LMC and 2% of the SMC; in these 
regions the BM C is complete for sources brighter 
than 200 mJy ( Seale et al.l 2014ll . This is a fac¬ 
tor of two greater than the completeness limits in 
low surface brightness regions. If our sample is 
distributed evenly across the Magellanic Clouds, 
then bright extended emission in the Magellanic 
Clouds may mask around five bona fide evolved 
stars with far-IR emission greater than 100 mJy. 
If we also take into account sources in the BMC 
which may potentially be evolved but were not in¬ 
cluded in our final list (as they are more likely to 
be YSOs) then this estimate increases to approxi- 
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Fig. 12.— SPIRE F 250 and F 350 CMD for the con¬ 
firmed evolved sources (blue) in the LMC (top) 
and SMC (bottom) overlaid on a Hess diagram of 
the HERITAGE Band-Matched Catalog sources. 
The grey circles represent the initial list of BMC 
candidate evolved stars and the red squares show 
the candidate evolved stars which have been con¬ 
firmed spectroscopically as YSOs. The dashed line 
marks the completeness limit of the BMC at low 
backgrounds. 
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Fig. 13.— IRAC Fs and PACS Fioo CMD. De¬ 
scription as in Figure [T^ 
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mately 10 sources. The Magellanic Clouds should 
therefore host between 18 and 23 post-main se¬ 
quence stars with far-lR fluxes greater than the 
BMC completeness limit, rather than the 13 we 
currently detect. The remaining 22 sources in our 
sample are faint and fall below this limit; here we 
do not have enough statistical information to pro¬ 
vide a quantitative value for the total population 
of faint evolved stars with similar fluxes. 

The identification of evolved stars in the BMC 
is biased towards brighter sources; these are 
prominent above the background and have an SED 
which is well sampled. To be included in our final 
list of probable evolved stars, each object must 
also be detected as a point-source in the mid-lR 
and be previously identified as an evolved star can¬ 
didate in the literature or via photometric colour 
cuts. As such we do not discover any new objects 
in our vetted sample. Faint sources which are 
variable and optically bright (e.g. OGLE variable 
stars), but have a low mid-IR flux will be excluded 
from our list as they are difficult to characterize. 
To complicate matters further, faint sources have 
larger positional uncertainties in the far-IR and 
hence have a higher possibility of a photometric 
mismatch. 

In Figures [12] and |T3| we compare our final 
sample of evolved stars to the entire HERITAGE 
BMC, the BMC candidate evolved stars which 
were identified as mismatch or probable YSOs and 
the BMC evolved star candidates confirmed spec¬ 
troscopically as YSOs. As noted earlier, the spec¬ 
troscopically confirmed sources are limited to the 
brightest stars in our sample due to signal-to-noise 
ratio constraints. In the SPIRE F 250 and F 350 
CMD there is a considerable amount of overlap 
between the different populations, and the sources 
are primarily separated based on their brightness 
at 250 p,m. Here we see a clear distinction between 
the bright spectroscopically confirmed YSOs and 
the evolved stars at log(F25o)~2.7, a similar cut 
is seen for the SMC despite its lower dust content 
which suggests a robust upper limit to the 250- 
micron flux of evolved stars. 

The majority of the sources that we conclu¬ 
sively identify fall well above the completeness 
limit of the Herschel BMC. Although we detect 
a moderate number of faint evolved sources below 
this line we have invariably missed others; it is dif¬ 
ficult to put a number on how may more of these 


faint sources might be evolved. The small spread 
in the ^ 250/^350 colour and the limited number 
of evolved stars detected in both bands prohibits 
the classification of sources using only their far-IR 
flux. If we are to clearly identify evolved sources 
additional mid-infrared data is required. FigurefTSl 
shows a CMD which combines mid- and far-IR 
photometry; here the various classes of objects 
tend to occupy different regions of CMD space 
although we are not able to completely separate 
the two populations. Furthermore, by combining 
far-IR data with 8 ^m photometry we effectively 
eliminate dust clumps and faint background galax¬ 
ies from the CMD making it easier to distinguish 
the evolved sources. 

5. Summary and conclusions 

In this paper, we present a catalog of thirty 
five evolved stars and stellar end products which 
have been detected with PACS and SPIRE as part 
of the Herschel HERITAGE survey of the Magel¬ 
lanic Clouds. Our source identification method is 
based on mid-IR spectral confirmation, the SED 
characteristics, careful examination of the multi¬ 
wavelength images and includes constraints on the 
luminosity and comparisons with YSO models. 

The sources detected in our catalog cover a 
wide range in luminosity and hence initial mass. 
We have identified thirteen low- to intermediate 
mass post main-sequence stars, which are either 
currently losing-mass at high rates or are evolved 
post-superwind objects. We find ten high mass 
stars, including four of the fifteen known B[e] stars 
in the Magellanic Clouds, and three extreme RSGs 
which are highly enshrouded by dust. Far-IR emis¬ 
sion is also detected from the ejecta of SN1987A 
and from the supernova remnant N132D, both in 
the LMC. We have also identified far-IR emission 
from ten possible evolved star candidates. 

This catalog lays the groundwork for future 
studies of the far-IR emission from evolved stars 
and stellar end products in the Magellanic Clouds. 
In a future paper we will provide in-depth treat¬ 
ment of the far-IR emission via radiative transfer 
modeling to determine if the emission originates 
from the circumstellar envelope, or if it is due to 
swept-up ISM dust. 
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